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SUMMARY 


As a result of a study conducted by the Federal Aviation Administration in 1977 
and a report published in 1983 by the National Research Council, low-altitude wind 
shear has become recognized as an infrequent but significant hazard to all aircraft 
during terminal-area operations. A sensor has been developed for measuring the time 
rate of change of the specific total energy of an airplane with respect to the air 
mass, and the sensor was tested and evaluated on a De Ha vi Hand DHC-6 Twin Otter 
airplane and also on the Transport Systems Research Vehicle (TSRV) B-737 airplane. 
Since total energy rate has been shown to be useful in wind-shear detection and in 
basic longitudinal stability and control studies, this sensor may be applicable to 
the wind-shear problem. 

This paper presents a control system (both with and without energy-rate feed- 
back) designed for the approach to landing of the TSRV B-737 through a severe-wind- 
shear and gust environment to evaluate the application of the total energy-rate 
sensor to wind-shear penetration during the landing phase of flight. A plant-system 
model is developed which incorporates wind-shear dynamics equations with the airplane 
equations of motion, thus allowing the control systems to be analyzed under various 
wind-shear conditions. The control systems are designed by using modern control 
theory and are formulated by using the proportional-integral-filter/command-generator 
tracking (PIF/CGT) control-law configuration integrated with optimal-output feedback. 
The feedback designs are analyzed with frequency-domain control-theory techniques 
for the design case with no wind and also for an off-design case with wind shear. 
Control-system performance is evaluated by using a complete nonlinear simulation of 
the airplane in combination with a severe-wind-shear and gust data package. This 
performance evaluation is concerned with control-system stabililty and regulation 
capability only - operational performance aspects, such as 100-ft decision height 
criteria and touchdown dispersion, are not considered. The analysis and simulation 
results indicate very similar stability and performance characteristics for the two 
designs. This outcome may be due to the availability and accuracy of the measure- 
ments that were used in this study. Energy-rate feedback may be beneficial on com- 
mercial transports that are not equipped with highly accurate measurement systems 
(e.g., inertial navigation systems) which are capable of providing all the necessary 
information or on commuter and general aviation aircraft, which have fewer available 
measurements for use in automatic control systems. 


INTRODUCTION 

Wind shear (or wind variability) is characterized by wind changes that depend on 
both distance and time. It can therefore be thought of as a local variation of wind 
velocity with respect to distance which occurs at some particular time. The magni- 
tude and relative direction of a wind shear are quantified by dividing the velocity 
difference at two points by the distance between them, and this quantity is usually 
referred to as the wind-shear gradient. 

Low-altitude wind shear results from thunderstorms and other convective clouds. 
It is recognized as an infrequent but significant hazard to all aircraft, from small 
general aviation airplanes to larger jet transports, during the take-off and the 
approach or landing phases of flight. A 1977 Federal Aviation Administration ( FAA) 



study (ref. 1) reported on aircraft accidents or incidents related to low-altitude 
wind shear occurring from 1964 to 1975. During this period, 25 cases involving 
transport airplanes were identified in which low-altitude wind shear was a possible 
cause of accidents during either the take-off (2 cases) or the approach or landing 
(23 cases). The most disastrous and well-known of these cases - in fact, the case 
that motivated the FAA study - was the 1975 landing accident at John F. Kennedy 
International Airport, which resulted in 112 fatalities. A later report under the 
auspices of the National Research Council (ref. 2) identified 27 cases related to 
low-altitude wind shear during the years 1964 to 1982. Four of the original 
twenty-five cases were eliminated as wind-shear-related accidents or incidents, but 
6 new cases were identified for the time period 1976 to 1982, including the 1982 
take-off accident at New Orleans International Airport which involved 153 fatalities. 

Reference 2 lists 19 recommendations for reducing the flight hazards posed by 
low-altitude wind shear. Portions of two of the six recommendations listed under the 
section entitled "Aircraft Performance and Operations" are pertinent to the research 
presented in this paper. The following excerpts are from these recommendations: 

Wind-Shear Effects on Flight Characteristics 


The FAA should sponsor analytical and simulator investigations to determine the 
wind shear penetration and recovery capabilities of transport aircraft, based on 
various onboard detection, guidance, and control systems. 

Guidance and Control Aids 


Onboard sensors and guidance aids should be evaluated in a systematic manner to 
determine their merits for future development and for possible retrofit in 
existing aircraft. These include flight director modifications, ground speed/ 
airspeed flight mangement systems, vertical-acceleration sensors, and energy- 
rate sensors. 

When an aircraft encounters a wind shear, changes occur in the wind-velocity 
components along each of the axes of motion of the airplane. Longitudinal wind 
shears run along the fuselage of the aircraft and are generally referred to as head- 
wind or tailwind shears depending on the direction. Vertical wind shears result from 
updrafts and downdrafts, which occur along the vertical axis of the aircraft. Longi- 
tudinal and vertical wind shears can add energy to or subtract energy from an air- 
craft because they affect the aircraft speed and altitude (i.e., its kinetic and 
potential energy, respectively) (ref. 2). These shears therefore have a major effect 
on the aircraft phugoid motion, which is characterized by constant total energy with 
a continual exchange between the kinetic and potential energy components. Further- 
more, certain wind-shear conditions can be shown to actually drive the phugoid mode 
to instability. Since total energy is the sum of kinetic and potential energy, its 
time rate of change with respect to the air mass reflects changes in altitude and 
airspeed due to wind-velocity changes. Thus, if total energy rate can be measured, 
information relating to wind shear will be available. Also, a measurement that pro- 
vides relatively fast information concerning wind shear, such as total energy rate, 
is particularly desirable over measurements that depend on the slow response of large 
airplanes . 

Total energy rate has been used previously in research applications for wind- 
shear detection (ref. 3) and also for a basic piloted longitudinal stability and 
control study (ref. 4). In addition, a relatively simple and inexpensive total 
energy-rate sensor has been developed for measuring the time rate of change of the 
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specific energy (energy per unit weight) of an aircraft with respect to the air mass 
(refs. 5 to 9). The sensor system, consisting of a pressure-measuring probe in com- 
bination with an in-line acoustic filter and an altitude-rate transducer, was previ- 
ously evaluated on both a De Havilland DHC-6 Twin Otter airplane (ref. 10) and on the 
Transport Systems Research Vehicle (TSRV) B-737 airplane at the NASA Langley Research 
Center (ref. 11). Both studies indicated a possible application of the sensor to the 
wind-shear problem and particularly to the automatic control of an airplane on final 
approach. Also, since this sensor is simple and relatively inexpensive, it is a good 
candidate to be retrofit on existing aircraft. 

The purpose of this study is to determine the results of incorporating the 
energy-rate sensor into an automatic longitudinal control system for the approach to 
landing of a transport airplane in a severe-wind-shear and gust environment. Includ- 
ing gusts (or turbulence) with wind shear contributes to flight hazards and increases 
the chances of an accident (ref. 2). A system model is developed which incorporates 
the signal from the total energy-rate sensor (including sensor dynamics) into the 
airplane longitudinal equations of motion, and a control system is designed and ana- 
lyzed both with and without the energy-rate measurement. In order to account for 
wind disturbances in the control-system design and analysis, wind-velocity equations 
representative of a linearly varying wind shear are incorporated into the system 
model, and wind-gust components are added as forcing functions. The control struc- 
ture selected for the design presented in this paper is the proportional-integral- 
filter/command-generator tracking (PIF/CGT) (ref. 12) control-law formulation com- 
bined with optimal-output feedback. 

The feedback designs are analyzed with frequency-domain control-theory tech- 
niques, whereas control-system performance is evaluated with time-history plots from 
a complete nonlinear simulation of the airplane combined with a severe-wind-shear and 
gust data package developed by SRI International. This data package includes a 
reconstruction of the hazardous wind-shear environment that occurred during the 
Kennedy airport landing accident, and this profile is used in the evaluation. This 
study is concerned with control-system stability and control regulation capability 
when subjected to a particularly hazardous wind-shear and gust environment. Opera- 
tional performance aspects, such as 100-ft decision height criteria and touchdown 
dispersion, are important but were not considered in this research. 

This paper presents a description of the PIF/CGT control-law structure with out- 
put feedback, the development of the design model, and a summary of the design, anal- 
ysis, and simulation results. It is shown that the control-system concepts presented 
herein result in acceptable designs from the standpoint of stability and control reg- 
ulation, with or without use of the total energy-rate sensor, and that inclusion of 
the sensor is not particularly advantageous. This outcome may, however, be due to 
the availability and accuracy of the measurements that were used in this study. 
Energy-rate feedback may be an advantage on commercial transports that are not 
equipped with highly accurate measurement systems, or on commuter and general avia- 
tion aircraft, which typically have fewer available measurements for use in automatic 
control. 


DESIGN MODEL DEVELOPMENT 

In order to evaluate the application of the total energy-rate sensor to the 
wind-shear problem, a longitudinal control law was developed for the approach to 
landing of an airplane through wind shear. The control structure selected for the 
design was the proportional-integral-filter/command-generator tracking (PIF/CGT) 
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discrete control-law formulation (ref. 12) combined with optimal-output feedback 
(ref. 13). This section provides a brief description of the control structure 
followed by the models required for the design. 

The PIF formulation is a linear feedback control structure that allows the feed- 
back of proportional and integrated states and that incorporates a control command 
filter to accommodate computational delay and to provide high-frequency-noise sup- 
pression. Proportional plus integral feedback provides partial trim accommodation in 
that the trim for the integrated states can be maintained while the proportional 
states go to new trims. The control command filter is implemented through the use of 
commanded control rates combined with proportional control-position/command feedback. 
The use of output feedback in the design precludes the necessity of estimating states 
so that only measurable aircraft parameters need be included in the output, and sen- 
sor, actuator, and engine dynamics can be accounted for in the state model. The PIF 
feedback design model is therefore comprised of a linear state (perturbation) model 
which incorporates plant, control, and integrator states in the state vector and uses 
control rates as commands. The output vector consists of selected measurable air- 
craft parameters in addition to the control and integrator states. The feedback 
gains are computed for the discrete system model in order to obtain a direct digital 
design. 

The CGT formulation (ref. 14) provides a means of incorporating command- 
generating and tracking capabilities into the control-system design through the use 
of feedforward control gains. The feedforward design is implemented independently of 
the feedback design and requires a system design model in addition to a linear com- 
mand model for the system to track. The system design model consists of the plant- 
state model used in the PIF design with the integrator states as its output, and the 
command model is formulated such that its output also corresponds to the integrator 
states. This construction enables the commands to these states to be tracked with 
type 1 system properties, thus enhancing the trim accommodation aspect of the design 
provided by the integrator states. 

The PIF/CGT feedback and feedforward structures are combined in an incremental 
form so that the flight equations can be expressed in terms of total quantities 
rather than perturbations. This precludes the necessity of storing trim values. A 
summary of the equations involved in the PIF formulation, in the optimal-output feed- 
back design, in the command generator for the feedforward design, and in the integra- 
tion of these structures for flight implementation are presented in appendix A. 

The TSRV B-737 airplane was selected for the design of this study for several 
reasons. First, computer programs were available in-house for determining the 
linear-perturbation model for some nominal flight condition (ref. 15) and also for 
performing a full nonlinear simulation of the airplane dynamic behavior in the pres- 
ence of sensor noise and external disturbances such as wind shear. Second, the 
energy-rate sensor system had already been installed, calibrated, and evaluated on 
this airplane in a previous study (ref. 11), which yielded a fairly accurate second- 
order dynamic model approximation to use in the design. Also, this airplane, with 
the sensors already installed, could be used for flight tests, if indicated. The PIF 
augmented feedback-design model, including the plant-, control-, and integrator-state 
models, is now presented, followed by the command model used in the feedforward 
design. 
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Feedback-Design Model 


The continuous plant-, control-, and integrator-state models required for the 
feedback design are presented first, followed by the PIF augmented design model. 

Plant-state model ,- The linearized continuous plant-state model combines the 
longitudinal aircraft states with those required for modeling the total energy-rate 
sensor (including sensor dynamics). Wind equations are also incorporated into the 
plant model so that wind disturbances can be accounted for in the control-system 
design and analysis. Other states included in this model represent altitude and 
engine dynamics. Appropriate output parameters were selected to be fed back to the 
controls. The plant-state model is summarized below, and all variables (unless 
otherwise noted) represent perturbations: 


X = AX + BD + B w W ( 1 ) 

Y x = H X X (2) 


where 


X = [u, w, q, 0, w w , h e f , h e>t , X th , h] T (3) 

(4) 

(5) 

Y X = t V a' 0 ' h ' Vtf {6) 


t fi e' *tJ 


W = [u , w l 1 
L w,g w,g 1 


The aircraft longitudinal states are given by u, w, q, and 0 in equation (3) • 
These represent, respectively, the longitudinal- and vertical-body-axis inertial 
velocity components, the pitch rate, and the pitch attitude. The longitudinal air- 
craft controls include elevator and throttle positions, and these are represented as 
6 e and 6^ in equation (4), The linearized aerodynamic coefficients, needed for 
the A and B matrices of equation (1), were obtained from the computer program 
described in reference 15 for the flight condition listed in table I, 

The next two states shown in equation (3) represent the longitudinal and verti- 
cal wind-velocity components u^^ and w w , expressed in the positive longitudinal and 
vertical directions of the body axes. These terms were included as states in the 
design model so the control law could be designed and analyzed for wind-shear dis- 
turbances. A linearized state wind model was therefore derived for incorporation 
with the airframe dynamic model. The wind model describes a wind-shear velocity for 
which the horizontal and vertical components vary linearly with respect to altitude. 
(Wind-velocity dependence on along-track position was not accounted for in this 
model, but it could easily be incorporated into the model for future studies.) Wind 
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velocities are obtained by integrating horizontal and vertical wind accelerations • 
Since the optimal-output feedback algorithm requires a strictly stable model, how- 
ever, the wind-model equations used in the design represent a low-pass-filter approx- 
imation to this integration with poles located slightly into the left half-plane. 

The poles were subsequently placed back at the origin to represent a true integration 
for the analysis. This wind-shear model was patterned after one that was presented 
in reference 16, and a complete derivation is provided in appendix B. 

Wind gusts are incorporated into the plant-state model as an additional forcing 
function, which is represented as W in equation (1). As defined in equation (5), 

W is a two-element vector containing the longitudinal and vertical wind-gust compo- 
nents. The wind-gust disturbance matrix B w relates these components to the air- 
craft states through the same relationships as the aircraft velocity components 
(ref. 17). This wind-gust vector was added to the design model so that wind gusts 
could be incorporated into the optimal-output feedback design. 

. • 

The next two states that appear in equation (3), h 0 ^ and h Q t , represent the 
second-order dynamic model of the energy-rate sensor. The model consists of a lin- 
earized expression for the specific total energy rate of the airplane h e passed 
through a second-order low-pass filter with appropriate time constants (ref. 11) to 
model the energy-rate sensor dynamics. Appendix C contains a complete derivation of 
the linearized energy-rate equation, given by equation (Cl 9), and the linear state 
model used for the sensor dynamics, given by equation (C20). 

The engine-dynamics state X t ^ was modeled as a first-order lag with a cutoff 
frequency set equal to 0.75 rad/sec, which is the numerical mean of the cutoff fre- 
quencies associated with the engine response during acceleration and deceleration. 
This value was used so the average engine response between its two operational modes 
would be reflected in the design. 

The altitude state h was modeled with the equation for altitude rate given in 
equation (Cl 5) of appendix C, since potential energy rate per unit weight is equiva- 
lent to altitude rate by definition. 

The output variables selected for the design are defined in equation (6). These 
include airspeed V a , pitch rate q, pitch attitude 9, altitude h, altitude rate 
h, and total energy rate h e t# Only measurable aircraft quantities were included as 
outputs in the design model . 9 The elements of the A, B, B w , and matrices of 

equations (1) and (2) are defined in appendix D. 

Control-state model .- The control-state model required for the control command 
filter is formulated for the control-rate equation by 


U = V 


(7) 


where 


V = 




(8) 
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Integrator-state model .- The states that were integrated in the PIF design model 
were altitude and airspeed. Altitude was integrated with the plant altitude state, 
and airspeed was integrated with the coefficients defined in equations (D9) to (D13) 
of appendix D. The integrator-state model can be written as 


Z = AX + B U 
z z 


(9) 


where 


Z = [/ h, / v a f 


( 10 ) 


and the matrix elements of A z and B z are defined in appendix D. 

Augmented PIF-state model .- The complete augmented-f eedback -design model was 
obtained by combining the plant-, control-, and integrator-state models developed in 
the preceding three sections. This model is summarized here as follows: 


• 





X 


A 

B 

0 

• 





o 

= 

0 

0 

0 

• 





z 


A 

B z 

0 


X 


0 


B 




w 

D 

+ 

I 

V + 

0 

Z 


0 


0 


(11 ) 



Y 


U 



0 



where 


T o - [««.. «thf 


( 12 ) 


( 13 ) 


*, -[/*•/ V ,1 T 


The terms H^j and H z of equation (12) are two-dimensional identity matrices, which 
enable the control-position commands and integrated states to be fed back to the con- 
trols. The above continuous state model was then discretized for a sampling period 
of 0.1 sec using the ORACLS computer program (ref. 18) as described in the n PIF 
Formulation” section of appendix A. 
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Command Model for Feedforward Control Design 

The linear discrete command model that was used for the feedforward control 
design contained one state and two outputs for commanding altitude and airspeed (the 
same aircraft states being integrated) and is given by 


< h «'k + I4T 01 


V 


a,c 


h 


~1 " 


0 

0“ 


h„ 

m 







c 


= 


[h m ] k + 





V 


_0 _ 


_0 

1 _ 


V 

a,m _ 

k 






a,c 


(1 5) 


(16) 


where h and V are the desired-altitude-rate and desired-airspeed inputs to 

C u f c 

the command model, and h m and V a are the commanded altitude and airspeed model 
outputs to be followed. The discrete sampling period AT is equal to 0.1 sec. 


CONTROL-SYSTEM DESIGN AND ANALYSIS 

The model described in the preceding section was used to design the PIF/CGT con- 
trol system for the approach to landing of the TSRV B-737 airplane. This section 
presents the feedback designs, with and without energy-rate feedback, that were ob- 
tained with the output feedback-design program (ref. 13), an analysis of each design 
from the computer programs of reference 19 and an in-house singular-value package, 
the feedforward design which was obtained in the manner discussed in reference 14, 
and an evaluation using a nonlinear simulation program of the TSRV B-737 airplane. 


Feedback Design and Analysis With Total Energy-Rate Feedback 

In order to compute the feedback gains via the optimal-output feedback -design 
program (ref. 13), an initial set of weights for the states and controls and an ini- 
tial set of process- and measurement-noise uncertainties had to be assigned so that 
the discrete weighting and covariance matrices required by the output feedback algo- 
rithm could be obtained. (See appendix A.) The state and control weights were 
obtained by computing the inverse square of the maximum allowable deviation for each 
state and control. These deviations were selected based on the corresponding trim 
values. The process-noise uncertainties were obtained by selecting reasonable values 
with respect to trim (based on engineering judgement) for an initial condition error 
for each of the states and controls and by assigning uncertainty values for severe 
wind gusts. Measurement-noise uncertainties were selected to represent typical noise 
spectra for the sensors that are onboard the TSRV B-737 airplane. 

The initial set of discrete weighting and covariance matrices was then input to 
the optimal-output feedback program along with the discrete system matrices (for the 
no-wind-shear case) in order to obtain an initial set of feedback gains. These gains 
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were then evaluated with frequency-domain analysis techniques and nonlinear simula- 
tion time histories. Weightings or uncertainties or both were modified to improve 
system stability and performance characteristics , and a new set of gains was computed 
and evaluated. This process continued until a final set of gains was obtained. A 
listing of all the weightings and uncertainties used in the final design is presented 
in table II, and the corresponding feedback gains are listed in table III. 

In order to analyze each design, the closed-loop eigenvalues were computed, and 
the gain and phase margins were obtained using Bode and singular-value plots. The 
Bode plots were computed with the loop broken at the actuator inputs (rather than at 
the rate commands) for each control separately while the other control loop remained 
closed. Only the input corresponding to the opened control response was considered. 
The singular-value plots were computed with both control loops broken at the actuator 
inputs for additive and multiplicative disturbances. This analysis was performed for 
the design case with no wind shear and also for a selected case with severe wind 
shear, both of which were run with the small wind poles (referred to previously) set 
to zero to represent the true wind-shear integration. The gradients for the off- 
design case were estimated from a reconstruction of the Kennedy airport wind-shear 
environment (referred to in the "Introduction"), which was developed as a software 
package by SRI International and is shown in figure 1 . This profile was generated 
for the nominal flight path of the airplane under control of the above design, with 
flare being initiated at 42 ft above the runway. The altitude and along- track posi- 
tion shown in the figure are plotted with respect to the airplane center of gravity 
(c.g.), and the zero point on the x cg scale is located at the end of the runway 
before touchdown. The gradients were calculated for specific time periods based on a 
nominal sink rate of 11 ft/sec. As shown in figure 1, the most severe gradients of 
the horizontal-wind-shear component occurred approximately from 46 to 54 sec and from 
74 to 83 sec, and those for the vertical-wind-shear component occurred from 56 to 
60 sec, from 64 to 74 sec, and from 74 to 86 sec. These gradients were then paired 
together to represent three realistic but severe environments. The open-loop eigen- 
values were computed for each, and the worst case was selected for the off-design 
analysis. Table IV lists the open-loop short-period and phugoid poles for each of 
these wind-shear environments, and table V summarizes the analysis results for the 
final-design and off-design cases. An increasing headwind (tailwind) is represented 
as +HW ( iTW) , and an increasing updraft (downdraft) is represented as fUD (+DD) . 

Since the phugoid poles are the most unstable for the third wind-shear condition 
listed in table IV, this was selected as the off-design worst case condition for the 
feedback analysis. The eigenvalues given in table V are listed with the phugoid 
poles first, the short-period poles second, followed by the remaining poles of the 
system listed in order of increasing value. In the open-loop eigenvalues for the 
design case, -1.13 and -2.40 represent respectively the acoustic-filter and trans- 
ducer dynamics of the energy-rate sensor, -0.75 represents the engine-dynamics state, 
two of the poles at the origin represent the wind states (since the small wind poles 
were set to zero to represent a true wind-shear integration for the analysis), and 
the five remaining poles at the origin represent the altitude state, the elevator 
control state, the throttle control state, and the two integrator states for altitude 
and airspeed. The unstable pole listed under the closed-loop eigenvalues for the 
off-design case is due to the wind states (which are uncontrollable) and is not 
indicative of an unstable aircraft system. The gain and phase margins (GM and PM) 
computed in the Bode analysis represent the stability margins of the system when one 
control loop is opened. The gain and phase margins computed in the singular value 
analysis, however, represent the guaranteed stability margins of the system for 
either additive or multiplicative disturbances when both control loops are open. 
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Feedback Design and Analysis Without Energy-Rate Feedback 


The feedback system presented in the preceding section was redesigned without 
energy-rate feedback with the same weights and uncertainties as those given in 
table II. The values associated with the energy-rate signal, however, were modified 
such that energy rate could be removed from the system without reducing the order of 
the model. These modifications included reducing the weighting on the two energy- 
rate states to zero, decreasing the uncertainty associated with these states to 
0.001, and increasing the uncertainty associated with the energy-rate measurement to 
100. All these changes tended to force the gains assigned to the energy-rate signal 
toward zero, which effectively removed it from the feedback loop. Table VI lists the 
gains obtained for the redesign without energy-rate feedback, and the corresponding 
analysis is summarized in table VII for both the design and the off-design case. 

The open-loop eigenvalues are not repeated in table VII because they are the 
same as those listed in table V. The analysis results indicate some differences in 
the closed-loop eigenvalues and stability margins for the two designs, but most of 
these differences are small. For the off-design case, feeding back energy rate 
appears to move the phugoid poles slightly farther into the left half-plane. In 
addition, the phase margins for the energy-rate-feedback design are slightly better, 
particularly when the throttle loop is open. These results suggest that the energy- 
rate signal may provide added phugoid-mode stability in wind shear, but the amount 
may be negligible in this application. Removing the energy-rate signal, alterna- 
tively, seems to provide a little more damping to the phugoid and short-period modes 
for both the design and the off-design case. 


Feedforward Design and Analysis 

The feedforward gains were computed as in reference 14 for the altitude and air- 
speed command model defined in equations (15) and (16). Since the state model and 
the states to be integrated were not altered for the two feedback designs (with and 
without energy-rate feedback), the feedforward design was the same for both systems. 
The feedforward gains are therefore summarized for either system as follows: 
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In order to analyze the feedforward control system, a nonlinear simulation of 
the TSRV B-737 airplane was used with implementation of altitude and airspeed com- 
mands which were computed with the command model of equations (15) and (16) for the 
altitude rate and airspeed input commands shown in figure 2. The altitude-rate and 
airspeed command-model inputs were computed from a constant desired-altitude acceler- 
ation and airspeed rate, respectively. Figure 3 shows the results of the nonlinear 
simulation of the TSRV B-737 following these commands in an environment which is 
uncorrupted by either noise or wind. Since no wind was introduced into the system, 
the ground-speed plot of figure 3 shows the airplane tracking the velocity (or air- 
speed) command. The altitude-rate response of the airplane is also following the 
desired-altitude rate. The plots of airspeed error AV a and altitude error Ah, 
which indicate the difference between the actual output and the commanded value, show 
that the airplane altitude and airspeed followed the command-model outputs except 
during the transitions, with airspeed having the slower response. The altitude mea- 
surement used in this simulation run was radar altitude. This signal, however, is 
generally replaced on the airplane by a third-order complementary-filter estimate 
(ref. 20), which is computed with respect to the runway coordinates. Figure 4 shows 
the same command-fol lowing simulation run with radar altitude replaced by this alti- 
tude estimate from a complementary filter. During flight, this filter is initialized 
prior to engaging the control system so that initialization transients can damp out. 
This takes approximately 60 to 80 sec. In this simulation, however, the complemen- 
tary filter was not initialized prior to engaging the control system. Therefore, in 
the figure, it is clear that the altitude response has been degraded due to these 
initialization transients introduced by the complementary filter, particularly since 
the filter dynamics were not included in the design model. This effect is evident 
during the first 60 sec of the altitude-error plot, which is consistent with the time 
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constants associated with the complementary- filter response. The simulations pre- 
sented in the next section all use this complementary-filter estimate for altitude 
rather than radar altitude. Therefore, additional error in the Ah plots will be 
present, particularly during the first 60 sec, that should not occur in actual 
practice . 


CONTROL-SYSTEM PERFORMANCE 

In order to evaluate control-system stability and control regulation capability 
under hazardous wind-shear and gust conditions, the control-system designs were 
implemented (using the flight equations presented in appendix A) in a nonlinear simu- 
lation of the TSRV B-737 which was used in conjunction with the SRI International 
wind-shear and gust data package. Operational performance aspects, such as 100-ft 
decision height criteria and touchdown dispersion, were not considered in this evalu- 
ation. The simulations began with the trimmed airplane located on a 3° glide slope 
at a mean-sea-level altitude of approximately 1300 ft. The airplane was commanded to 
descend down the glide slope at constant airspeed until the airplane initiated flare 
at 42 ft above the runway. During flare, the aircraft was commanded to follow an 
exponential path with a constant deceleration of 3 ft/sec 2 . The simulations ended 
when touchdown occurred. The airplane equations of motion were computed every 
0.05 sec, and the control commands were updated every 0.10 sec for a desired altitude 
along the glide slope and a desired airspeed. A slower sampling rate was sufficient 
for implementing this control system because it is a direct digital-control-law 
design. This is a major advantage of performing a direct digital design over dis- 
cretizing a continuous control-law design. The 0.10-sec time increment was selected 
for this application based on previous experience in implementing a digital control 
system on the TSRV B-737 (ref. 21). All the simulation results discussed in this 
section were for the Kennedy airport wind-shear reconstruction shown in figure 1 
combined with a corresponding gust package that contains root-mean-square (rms) gust 
intensities up to 1 3 ft/sec (ref. 22). The lateral-wind-velocity component was set 
to zero, however, to reduce lateral effects. Figure 5 shows a time-history plot of 
this wind-shear and gust environment as it was used in the simulations, although time 
varies slightly for each run due to varying ground speed. All the design runs pre- 
sented used the altitude estimate from the third-order complementary filter rather 
than the radar altitude, and measurement noise which was representative of the sen- 
sors onboard the airplane was incorporated into the simulations. Position and rate 
limits for the control-surface servos and actuators were also included. Lateral con- 
trol was implemented within the simulation for maintaining the wings level and the 
desired heading, but there was no closed-loop control for localizer tracking. This 
section presents the simulation results for the control-system design both with and 
without energy-rate feedback (figs. 6 to 9), followed by the performance of the TSRV 
B-737 instrument landing system (ILS) autoland system in this environment as a com- 
parison (fig. 10). In order to provide information regarding aircraft altitude in 
these results (figs. 6 to 10), the point at which the airplane is approximately 
100 ft above the runway is indicated along the time scale of the figures. A brief 
discussion of robustness is then presented for several off-nominal cases in which 
basic airplane parameters such as c.g. location, weight, and airspeed were varied. 
Slight variations in these parameters had little effect on control-system perfor- 
mance, so inclusion of these effects as separate figures was not warranted. However, 
in order to present these effects in this paper, the results presented in figures 6 
to 9 are those for which the airplane c.g. was slightly different from the design. 

The glide-slope angle was then increased to 5° and these results for each design are 
presented in figures 11 and 12. All these plots represent total aircraft quantities 
except where noted. 
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Performance of the Design With Energy-Rate Feedback 

Nonlinear simulation results for the control-system design with total energy- 
rate feedback are shown in figures 6 and 7. The first case (fig. 6) contains air- 
speed in the measurement set (see eq. (6)), whereas the second case (fig. 7) contains 
the velocity V, which distributes the velocity feedback gains between airspeed 
and ground speed V g according to the following implementation: 

v k ■ <’ - °> v a,k * ° v g,k < 21) 


where 


c = 0 


(V a,» < V 


( 22 ) 


is for total airspeed control 


c 




< V 


a,m 


< V + V) 

g l 


(23) 


is for linear distribution between airspeed and ground speed, and 


c = 1 


(V 


a,m 


v g + 




(24) 


is for total ground-speed control. 


In equations (22) to (24), V a m is the airspeed generated by the CGT model, and 
is determined by the allowable ' lower limit on ground speed. This modification does 
not alter the validity of the design analyses since the total velocity gains remain 
unchanged. in most cases, V a will have the larger component of feedback gain, and 
for the case of a tailwind, all the feedback gain is applied to regulation of V a . 
The limit of 50 ft/sec was selected for to prevent stall in case of a sudden 

wind change from a headwind to a tailwind. 

From figure 6, the worst altitude error Ah is -65 ft, with one- third of this 
error caused by the complementary-filter initialization transients. The airspeed 
data show tight tracking through the peak headwind of approximately 60 ft/sec, which 
occurred at around 70 sec into the run. During the time period of 85 to 1 1 0 sec, 
however, AV decreases by 70 ft/sec, V drops to 130 ft/sec, a increases to 
17°, 0 increases to 24°, and 6 th reaches saturation (60°). It is interesting to 

note that this time period corresponds to a decreasing headwind and a decreasing 
downdraft, as shown in figure 5, and can be related to the increasing tailwind and 
increasing updraft condition evaluated in the "Control-System Design and Analysis" 
section and found to induce phugoid-mode instability. Although the simulation does 
not indicate that stall occurred, the aircraft may actually reach stall under these 
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conditions due to several factors that were not accounted for in the simulation. For 
example, slight asymmetries that exist in the wings of the aircraft were not modeled 
in the simulation. Also, lateral effects due to winds and closed-loop localizer 
tracking were not included. Either of these factors could cause a wing to stall 
under the conditions of low V a and high a shown in figure 6, resulting in the 
aircraft reaching stall. 

The implementation of the control system was therefore modified as previously 
discussed to limit decreases in V g while maintaining enough V a control to prevent 
stall during severe shifts from headwind to tailwind shears. These results are pre- 
sented in figure 7. In contrast to figure 6, this figure shows that V a is allowed 
to increase during the large headwind, with the subsequent decrease being much 
smaller than in the previous case. In fact, never decreases very much below the 

commanded value. Furthermore, just touches saturation, 0 remains relatively 

low during the severe-wind-shear time period, and a remains relatively constant 
throughout the entire approach and landing. The total run time is approximately 
5 sec shorter because of better regulation of V g . 

Operational aspects that were not addressed in this study include the 100-ft 
decision height criteria for continuing the approach or implementing a go-around 
(this includes maintaining an adequate thrust reserve for effecting a go-around, if 
necessary) and touchdown conditions for an acceptable landing. These operational 
aspects could probably be improved by using additional implementation techniques 
similar to the one presented here and by gain scheduling. For example, an upper V g 
limit was implemented in the same manner as the lower limit in an effort to reduce 
the high V g at touchdown, and this proved to be effective in accomplishing the 
desired reduction. These results were not included in the figures, however, because 
they were preliminary in nature and additional modifications would still be required 
for an operationally acceptable touchdown, which was not an objective of this study. 


Performance of the Design Without Energy-Rate Feedback 

The same simulation results as those just presented for the design with energy- 
rate feedback are repeated in figures 8 and 9 for the redesign of the control system 
without energy-rate feedback. As shown in figure 8, removing the energy-rate mea- 
surement actually appears to slightly improve control regulation capability. The 
value AV a is reduced to approximately -56 ft/sec, a is reduced to 15°, and 0 is 
reduced to about 22°, while the maximum Ah error remains about the same. In addi- 
tion, and activity is slightly reduced when the energy-rate signal is 

removed. The results shown in figure 9 also indicate a slight improvement in the 
control-regulation capability of the redesign without energy-rate feedback, although 
the differences between this figure and figure 7 are less pronounced than those 
between figure 8 and figure 6. The overall results for the two designs are very 
similar, however, and both sets of results indicate a dramatic improvement in the 
regulation capability of the implementations with distributed V a and V g control 
over those with V a control only. 


TSRV B-737 ILS Autoland-System Performance 

The performance of the instrument landing system (ILS) autoland control system 
onboard the TSRV B-737 airplane (ref. 20) is presented in figure 10 as a baseline run 
for comparison. The initial conditions and trim values for this run are set equal to 
those which were used for the previous design runs. The ILS system is an analog 
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control law that is implemented digitally on the airplane flight-control computers at 
a rate of 20 times per second, and the simulation program has an equivalent implemen- 
tation. Figure 10 shows that the ILS system is able to maintain acceptable stability 
and control regulation until the airplane encounters the decreasing-headwind/ 
decreasing-downdraf t condition mentioned previously. During this time period, AV a 
decreases by 65 ft/sec, Ah decreases by 138 ft, a increases to 24°, 0 increases 
to 30°, and 6 e reaches its lower saturation limit. The value of during this 

period touches 60°. The 24° limitation on a shown in figure 10 is imposed within 
the simulation because of the unavailability of aerodynamic data beyond this point. 

In actuality, a would continue to increase until the airplane reaches stall. The 
simulation indicates that the airplane may begin to stall during this period, but 
that altitude is recovered (and, in fact, reaches 256 ft above the glide slope) sub- 
sequent to 5 e saturating its upper limit. A definite stall condition probably 
would occur in reality, however, without the limitation on a. Also, wing asymme- 
tries and lateral wind effects could also combine with this condition, as discussed 
previously, to cause stall during this period. Furthermore, if the decreasing- 
headwind/decreasing-downdraf t condition had occurred any later in the run, the air- 
plane probably would have impacted the ground even though a recovery may have been 
theoretically possible. The airplane lands in the simulation, however, and the large 
positive Ah error indicated at touchdown results from the glide slope extending 
beneath the runway beyond the glide-slope intercept point. 

A comparison of the ILS results of figure 10 with those shown in figures 7 and 9 
indicates that both designs are able to maintain much better control stability and 
regulation quality in this environment than the ILS system. During the same critical 
time period, the worst error occurs in 0, which increases to 11° and 10° respec- 
tively for figures 7 and 9, and just touches its 60° position limit. All other 

parameters, however, are held fairly constant throughout this period, and stays 

below 6.5° in figure 7 (design with energy-rate feedback) and below 4.7° in figure 9 
(redesign without energy-rate feedback). Furthermore, even the design results pre- 
sented in figures 6 and 8 for the initial implementation are preferable to those 
shown in figure 10 for the ILS system. This may be because of the optimal integrated 
elevator-throttle structure of these designs, combined with the ability to account 
for process noise, measurement noise, and control-position and rate saturation limits 
during the design process. 


Robustness 

A precursory evaluation of the control-system designs for robustness was ob- 
tained for three off-design cases by varying the c.g. location, weight, and airspeed 
of the aircraft and repeating the simulation runs with the previously described 
velocity gain distribution for each design. The off-design cases selected for this 
analysis are presented in the following table: 


Case 

c.g., percent MAC 

Weight, lb 

Airspeed, knots 

1 

10 

85 000 

135 

2 

20 

85 000 

125 

3 

30 

75 000 

119 


15 




The effect of this variation on the simulation results was not significant enough to 
warrant presenting these runs as separate figures. However, in order to include the 
effect of slight parameter variations on the performance of the control-system 
designs presented in this paper, the preceding simulation results shown in figures 6 
to 9 are for off-design case 2 listed above, which represents a slightly off-nominal 
c.g. location. (The aircraft parameters for the design included a c.g. location at 
25 percent MAC rather than at 20 percent MAC.) This off-design case was then re- 
peated for each design with the glide-slope angle increased from 3° to 5°, and these 
results are presented in figures 11 and 12. It is shown in these figures that the 
pitch attitude and throttle activity are worse in these runs, but the general perfor- 
mance of both designs is very similar to the results presented in figures 7 and 9. 


SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 

A control system with and without energy-rate feedback has been developed for 
the approach to landing of the Transport Systems Research Vehicle (TSRV) B-737 air- 
plane through a severe-wind-shear and gust environment to evaluate the application of 
a total energy-rate sensor to the wind-shear problem. The analysis and simulation 
results for the two designs can be summarized as follows: 

1 . A simple implementation technique for distributing the velocity feedback 
gains between airspeed and ground speed proved to be very effective in reducing pitch 
attitude, angle of attack, and throttle activity during severe-wind-shear and gust 
penetration while maintaining good airspeed and altitude control. This implementa- 
tion technique does not alter the design or the validity of the analysis since the 
total feedback gains remain unchanged, but operational performance criteria may be 
affected. 

2. The control-law design, both with and without total energy-rate feedback, 
demonstrated the capability of following both altitude and airspeed commands with 
zero steady-state error and of penetrating a severe-wind-shear and gust environment 
during the approach to landing. 

3. Integrating the airplane and wind-shear dynamics as coupled differential 
equations indicated wind-shear conditions for which the airplane phugoid mode was 
unstable. This modeling approach allowed the control law to be analyzed for these 
unstable conditions. 

4. Minor differences appeared in the analysis and simulation results of the 
control-law design with total energy-rate feedback compared with those of the 
control-law design without total energy-rate feedback, but none of these were of 
great significance. 

5. The instrument landing system (ILS) autoland system onboard the TSRV B-737 
airplane was much less effective in penetrating the severe-wind-shear and gust 
environment than were the control-law designs (with and without total energy-rate 
feedback) considered in this paper. 

From these results, the following conclusions may be drawn: 

1 . The energy-rate sensor may be better applied to aircraft with a more limited 
and less accurate measurement field, such as for commercial transports that are not 
equipped with an inertial navigation system or for commuter and general aviation 
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aircraft, since the other measurements studied herein appear to provide all the 
necessary information. 

2. The design concepts, control-law structure, and implementation techniques 
presented in this paper are applicable, with or without use of the energy-rate 
sensor, to the wind-shear problem, and these designs are much more effective in main- 
taining stability and control regulation during severe-wind-shear and gust penetra- 
tion than conventional systems. This may be because of the optimal integrated multi- 
variable structure of these designs and the ability to account for process noise, 
measurement noise, and control-position and rate saturation limits during the design 
process. 

3. Incorporating wind-shear-dynamics equations with the airplane equations of 
motion for the plant design model is an effective means of providing control-system 
analysis capability for various wind-shear conditions. 

Based on the results and conclusions of this study, the following recommenda- 
tions are offered for future study: 

1 . A possible extension to this study would be to investigate the sensor for a 
reduced measurement field, such as for transports without the inertial navigation 
system or for a commuter or general aviation aircraft. The value of fast-response 
wind-related information provided by the total energy-rate sensor may be of greater 
importance there than in the application of this study. 

2. The wind-shear model presented in this paper includes wind component acceler- 
ation as a function of altitude rate. This wind-shear model should be expanded to 
include wind component acceleration as a function of the longitudinal velocity of the 
airplane. The expanded model would allow a more thorough investigation of sensitive 
wind-shear conditions and control-system analysis under these wind-shear conditions. 

3. Additional research should include operational design criteria, that is, 
maximum and minimum airspeed and ground-speed criteria, 100-ft decision height flight 
characteristics, and touchdown dispersion under various wind-shear and gust 
environments. 

4. A further application of these concepts might be to develop a similar design 
for take-off. 


NASA Langley Research Center 
Hampton, VA 23665 
January 18, 1985 
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APPENDIX A 


CONTROL-LAW STRUCTURE 

This appendix presents a summary of the PIF formulation, the output feedback 
equations, the CGT formulation, and the flight equations used in the design and eval- 
uation of the control systems presented in this paper. All variables represent 
perturbation quantities unless otherwise noted. 


PIF Formulation 

The continuous-time design model using the PIF structure is given by the 
following: 
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+ B W 

(A1 ) 
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where 


X 

u 

z 

V 

w 

Y 

Y x 

Y u 

Y z 

N 

A,B 

B w 

H 

A z' B z 


state vector for the plant 

control position-command state vector 

integrator-state vector 

control rate-command vector 

wind-gust vector 

output vector 

output vector for the plant states 
output vector for the control states 
output vector for the integrated states 
PIF output-noise vector 
plant-state and control-state matrices 
wind-gust disturbance matrix 
observation matrix 

transmission matrices for the integrator states 


The discrete form of the state equation given in equation ( A1 ) can be obtained 
by using the ORACLS computer program (ref. 18). The discrete design model can be 
written as follows: 


^ + 5v k * 5 »\ <fl8) 

\ ■ “k + \ <i9> 


where F, G, and G are the discrete matrices corresponding to A, B, and B 

w — - w 

as obtained from ORACLS. However, several approximations to F and G are made m 

order to implement both a zero-order hold for the control-rate commands and an Euler 

discrete integration (ref. 12). The final form for F and G, then, becomes 
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(A10) 
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0 


G = 


AT I 
0 


(A1 1 ) 


where AT is the discrete sampling period and I represents appropriately dimen- 
sioned identity matrices. 


Output Feedback Equations 

A stochastic, discrete optimal-output feedback design program (ref. 12) was used 
to compute the feedback gains for the selected outputs of equation (A9) . The 
continuous-time cost function J is given by 


(A12) 


Q is a diagonal state-weighting matrix with 

components given oy 


J = *{fo 


-T — T- l 

X QX + V RyV dt 


where E{ } is the expected value and 


Q = diag[Q x , Q 0 , Q z ] 


(A13) 


where Q x is the plant-state-weighting matrix, is the control-position command- 

state-weighting matrix, and Q z is the integrator-state-weighting matrix. The 
control-rate comma nd-weigh ting matrix is represented by the diagonal matrix R^. All 
weightings described in the "Control-System Design and Analysis" section are inserted 
into equation (A12). 

The equivalent discrete performance index, which is required in reference 13, is 
defined as 


J 



-T*- 


— T A 


T a 

V v ki 


(A1 4) 


A A — — 

where Q and R are the discrete weighting matrices corresponding to Q and R^ 

as calculated with the sampling period AT using reference 18. The term M is the 
discrete cross-weighting matrix between states and controls. 

The output feedback formulation also includes both plant process noise and mea- 
surement noise. Terms included in the plant process noise are randomly distributed 
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initial-condition errors, control-input random pseudonoise, and random wind-gust dis- 
turbances* The total discrete plant-process-noise covariance w is calculated as: 


w 


= FX.F T + f e^ t (BV B T + B W B T )e^ t 

1 J o L 1 P w g w' 


at 
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(A15) 

(A1 6a) 

( A1 6b) 
( A1 7a) 
(A1 7b) 
(A18a) 
(A18b) 


is the initial-condition error covariance, V is the control-input 


Analysis" section of this paper, 
to the discrete form using reference 18. 


where X. — 

1 _ P 

process-noise covariance, is the wind-gust process-noise covariance, and 

is the Kronecker delta function. The terms x^, v, and w contain respectively the 
state, control, and wind uncertainties described in the "Control System Design and 

The integral term in equation (A15) is transformed 
The measurement-noise covariance v is 

calculated in a similar manner as follows: 


e{n} = 0 


EflMT } = V 6. . 
1 1 m X] 


(A19a) 
(A1 9b) 


where N contains the measurement-noise uncertainties* 

The optimal-feedback gains K are directly related to the output measurement as 


V, = KY, 
k k 


(A20) 
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with the corresponding closed-loop transition matrix given by 

F = F + GKH ( A2 1 ) 

cl 


CGT Formulation 

The objective of command-generator tracking (CGT) is to cause selected outputs 
Y v of the airplane to optimally track the output Y _ v of a linearized command 
model by using feedforward control gains. The open-loop plant model can be written 
as 


x k+1 - + GO k 
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and the linearized command is given by 


^jk+l F m x m,k + G m°m r k+1 


(A24) 


Y m,k ^^jk + D m u m,k+1 


(A25) 


It should be noted that the output vector Y z ^ in equation (A23) corresponds to the 
airplane states being integrated in equation (A8) , and that the model outputs Vk 
in equation (A25) are the commanded values for these same states. The reason for 
using an advanced time step on the control is presented in reference 12. The feed- 
forward gains relate the model states ^ and model controls to ideal star 

trajectory X ,U as follows: 
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where each 


A i? 


represents a constant feedforward gain matrix. 


feedforward gams involves solving equation (A26) as follows: 
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A method for solving equation (A27) and the derivations involved are presented in 
reference 14. (A more detailed derivation is presented in TIM No. 612-2 of The 
Analytic Sciences Corporation entitled "Command Generator Tracking - The Discrete 
Time Case"; this document is not widely available.) 


Flight Equations 

The CGT model is integrated into the feedback design by letting the perturbation 
model (eq. (A8) ) represent the error between the plant perturbation states and con- 
trols and the star trajectory perturbation states and controls. If we use an incre- 
mental control-law form and rigorously substitute total quantities into the perturba- 
tion model (similar to ref. 12, which uses full-state feedback as opposed to the 
output feedback used in this paper) , the flight equations become 
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T T 

where p . p and P represent the rows of the optimal-output feedback solu- 
XZ T UZ ' ZZ **' 

tion matrix of equation (A14) that correspond to the integrator states. The sub- 
script T's in the above equations indicate these are total quantities. 
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APPENDIX B 


WIND-MODEL DERIVATION 

This appendix presents a complete derivation of the linear wind-shear states 
which were incorporated into the control-law design model presented in this paper. 

The wind-velocity vector in the inertial Earth-fixed reference frame is given by 


v „ = u 
w,E w,E 


W z 
w, E 


( B1 ) 


where x and z are unit vectors. This vector is then rotated to the body axes by 
using the following transformation: 
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The resulting body-axis expression for wind velocity is 


V = (U _ cos 0 - W „ sin 0)x, 
w w,E w,E 


+ (u w,E sin 6 + W w,E cos 0)z b 


(B3) 


Differentiating the and components of equation (B3) yields 


U w = U w,E COS 0 - U w,E 0 Sin 0 - W w,E sin 0 - W w,E 0 cos 0 


(B4) 


. . 
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for pitch rate in equations (B4) and (B5) yields 


U w = U w,E cos 0 - 2 U w,E sin 0 ' W w,E sin 0 - ^w,E COS 0 


" U w,E cos 0 " W w,E Sin 0 " W w 0 


(B7) 


W w = U w,E Sin 0 + 2 U w r E COS 0 + W w,E COS 0 
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= U WfE sin 0 + W WfE cos 0 + U W Q 
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The following expressions for U w E and W w E can now be obtained for a hypothet- 
ical linear wind shear that varies with respect to decreasing altitude: 
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some nominal constant horizontal wind velocity 
some nominal constant vertical wind velocity 
horizontal wind-velocity perturbation due to wind shear 
vertical wind -velocity perturbation due to wind shear 
horizontal wind-shear gradient, (ft/sec) /ft 
vertical wind-shear gradient, (ft/sec)/ft 
aircraft altitude, positive down (-H) 


Differentiating equations (B9) and (BIO) yields 
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where u „ and w _ represent the total wind perturbations from the nominal wind 
, . . W, E W, E 

velocities u _ and w tt . 

w,E,o w,E,o 
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. Equations (B1 1 ) and (B12) can be related to the airplane dynamics by defining 
Z E to be the sink rate of the airplane along the glide slope (i.e., for a frozen 
wind field) • The sink-rate equation as a function of the airplane longitudinal- and 
vertical-body-axis velocity components is given by 


-U sin 0 + W cos 0 


(B1 3) 


Substituting equation (B13) into equations (B1 1 ) and (B12) yields 


U w E = -U Z (U sin 0 - W cos 0) 


W w E = ~ w 2 ( u sin 0 - W cos 0) 


(B1 4) 


(B15) 


The complete nonlinear expressions for the longitudinal- and vertical-body-axis wind- 
acceleration components can now be determined by substituting equations (B14) and 
(B15) into equations (B7) and (B8) : 


U w = 


-U Z (U sin 0 - W cos 0) cos 


0 + 


W 2 (U sin 0 - W cos 0) sin 0 - W w Q 


(B16) 


W. = - 


U Z (U sin 0 - W cos 0) sin 0 - W Z (U sin 0 - W cos 0) cos 0 + U W Q 


(B17) 


The linearized expressions for the wind-acceleration components can be obtained 
by performing a perturbation expansion on equations (B16) and (B17) with the follow- 
ing substitutions: 


U = U Q + u 


(B1 8) 

> 

+ 

o 

II 


(B19) 

0 = 0 O + 0 


( B20 ) 

Q = q Q + q 


(B21 ) 

II 

O 

+ 

w w 

(B22) 

U w = ^,0 + 

% 

( B23 ) 
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where u Q , w Q , 0 O , and q Q are the nominal trimmed aircraft parameters, and u r 
w, 0, and q represent perturbations. Substitutinq equations (B18) to (B23) into 
equations (B16) and (B17), inserting appropriate trigonometric identities, expanding, 
and collecting terms yields the final form of the wind-model equations: 


K = [“ u 2 (0 * 5 sin 20 O ) + w z( sin2e o)] u + [ u z ( cos2 e o ) - W z (0.5 sin 20 o )]w + (-w^ 0 )q 

+ [— ( U u + W w ) cos 29 - (U w - W u ) sin 20 10 + (-q )w 

L z o zo o zo zo o J o w 

w w = [-U z (sin 2 0 o ) - W z (0.5 sin 20 Q ) ]u + [u z (0.5 sin 20 Q ) 

+ [~ (U z u o + W z w o ) sin 26 o + (U z w o - "zV cos 20 o] 6 + (< 3o )u w 

u w,o " “ u z [ u o (0 * 5 sin 20 o } ~ w o( cos20 o)] ‘ w z[ w o (0 ’ 5 sin 29 o ) - u o( sin20 o)] 

Vw , o 

w w,o = u z [ w o (0 * 5 sin 20 o ) - u o( sin20 o)] - w z[ u o (0 - 5 sin 20 o ) - w o( cos20 o)] 


5 sin 20 q ) ]w + 

( — w )q 

w, o ^ 

(-q )w 

( B24 ) 

o w 


W z( cos20 o)l w + 

(u w,o )c t 

(C to )U W 

( B25 ) 


(B26) 


+ *o u w,o 


( B27 ) 


Equations (B24) and (B25) represent the wind-shear perturbations from the desired 
flight path, and equations (B26) and (B27) represent the nominal wind-shear perturba- 
tions from the constant wind velocities ^\f f0 and along the desired flight 

path. The terms u^^ Q and w w Q are obtained by rotating the nominal constant 
wind-velocity components in the Earth-fixed inertial frame E Q an< 3 w w E o to 
their equivalent body-axis components with equation (B2) . The sign conventions 
implicit in the above equations are summarized as follows: 

u w E o > 0 for constant tailwind 
w w E o > 0 for constant downdraft 
U z > 0 for increasing tailwind with descent 
W z > 0 for increasing downdraft with descent 

Equations (B24) and (B25) are the linearized wind- state-model equations, and 
each wind-velocity component is obtainable through a pure integration of these accel- 
eration components. The wind-model equations used in the design, however, represent 
a low-pass-filter approximation to this integration (because of a restriction in the 
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control-law algorithm requiring a strictly stable system), and these equations are 
given below: 


u 


w 


[-U z (0.5 sin 20 o ) + W z (sin 2 e o )]u + [u z (cos 2 0 o ) 


+ ( — w )q + f-(U u + W w ) cos 20 - (U w - 

w,o L z o zo o zo 


+ ( -u ) u -f ( — q ) w 

w,w w o w 


- W z (0.5 sin 20 q ) ]w 

W u ) sin 20 10 
z o o J 

( B28 ) 


w 


[-U z (sin 2 0 o ) - W z (0.5 sin 20 Q ) ]u + [u z (0.5 sin 20 Q ) 

+ (u )q + f— (U u + W w ) sin 20 + (U w - W u ) 

w f o 1 zo zo o zo zo 


+ W z (cos 2 0 o )]w 


cos 20 0 

o J 


+ (q )u + ( — w )w 

o w w,co w 


( B29 ) 


where 


u o 

nominal inertial velocity 
X-direction of the body 

w o 

nominal inertial velocity 
Z-direction of the body 

^7,0 

nominal component of wind 
tion of the body axes 

w w,o 

nominal component of wind 
the body axes 


component of the airplane in the positive 
axes 

component of the airplane in the positive 
axes 

velocity in the positive longitudinal direc- 


velocity in the positive vertical direction of 


U z horizontal wind-shear gradient with respect to decreasing altitude 

W z vertical wind-shear gradient with respect to decreasing altitude 

Q nominal wind-shear acceleration in the positive X-direction of the body 
axes 

w w o nominal wind-shear acceleration in the positive Z-direction of the body 
axes 


u w,o) small low-pass cutoff frequency for the horizontal-body-axis wind- 
velocity component 

w w,w small low- pass cutoff frequency for the vertical-body-axis wind-velocity 
component 


The cutoff frequencies and w w ^ were selected to be small with respect to 

wind spectra and the system bandwidth* so that the frequency response of the wind 
"filters" would approximate an integrator for those frequencies. This approximation 
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was incorporated into the design model, because the output feedback algorithm used in 
the design requires that the initial system be strictly stable. Hence, incorporation 
of the small low-pass-filter poles adequately approximated the integrations while 
allowing the eigenvalues associated with the wind states to be located slightly into 
the left half-plane rather than at the origin, thus meeting the stability require- 
ment. These poles were subsequently set back at the origin for the control-system 
analysis . 
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APPENDIX C 


TOTAL ENERGY-RATE MODEL DERIVATION 

This appendix presents a complete derivation for the linear energy-rate equation 
used to represent the energy-rate sensor model presented in this paper. Total energy 
Erp can be expressed in terms of potential energy Ep and kinetic energy E^ as 
follows: 


E t = Ep + Efc = mgH + (1/2)mV^ 


(Cl ) 


where H indicates total altitude. The total specific energy per unit weight, or 
energy height H e , is obtained as 


H e = H + <1/2g)V* 


(C2 ) 


Differentiating equation (C2) yields the basic total energy-rate model used in the 
design: 


H e = H + (i/g)v a v a 


(C3) 


where the first term represents the potential energy-rate component of the airplane 
and the second term represents the kinetic energy-rate component. Since the poten- 
tial energy-rate component is simply the altitude rate, the body-axis equation can be 
expressed as 


Hp e = U sin 0 - W cos 0 


(C4 ) 


Since the kinetic energy-rate component was found in reference 11 to measure the air- 
speed and airspeed-rate terms along the longitudinal direction of the body axes, the 
kinetic energy-rate component is written as 


H ke = (1 /9 )V 


a,x 


a,x 


(C5) 


where 


V 


a, x 


U 


(C6 ) 


V 


a,x 


U 


w 


( C7 ) 
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Equation (C5) includes only the longitudinal-direction velocity and acceleration 
components because it was previously determined in reference 11 that this quantity 
provides a better representation of the total energy-rate sensor measurement. The 
term U in equation (C7) represents the airplane longitudinal acceleration and is a 
function of the airplane states and wind-velocity components in the longitudinal and 
vertical directions of the body axes. 

Equations (C4) and (C5) can be linearized through a perturbation expansion with 
the following substitutions: 


U = U Q 

+ u 

(C8) 

II 

£ 

+ w 

( C9 ) 

CD 

II 

CD 

0 

+ 0 

(CIO) 

V a,x 

V + v 

a,x,o a,x 

(C11 ) 

^a,x ~~ 

. . 

V + v 

a,x,o a,x 

(Cl 2 ) 

v a,x 

u - % 

(Cl 3 ) 

v a,x 

• . 

u - u 
w 

(Cl 4) 


Substituting equations (C8) to (CIO) into equation (C4) yields the linearized poten- 
tial energy-rate perturbation 


hp e = u sin 0 Q - w cos 0 Q + (u Q cos 0 Q + w Q sin 0 Q )0 


(Cl 5 ) 


and substituting equations (Cl 1 ) to (C14) into equation (C5) yields the linearized 
kinetic energy-rate perturbation 

\e = ( V a ,x,o/9 )u - ( ^a,x,o/9 )u w + (V a,x,o/^ “ (V a,x,o^^w (C16) 

where 


0 

K 

1 

(Cl 7 ) 

u w,o 

(Cl 8 ) 
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u = a n u + a 12 w + a 13 q + a 14 0 + + a 16 w w + b , , 6 e + b 12 6 th (C19) 

u w = w n u + w 12 w + w 13 q + w 14 0 + w 15Uw + w 16 w w + ^ f0 (C20 ) 


Adding equations (C15) and (C16), collecting terms, and substituting equations (C19) 
and ( C20 ) for the acceleration terms u and yields the final form of the 

energy-rate equation required for the design presented in this paper: 


h e = t sin 0 o + (v a,x,o/9) + ( v a ,x,o/9 )(a 11 " W 1 1 > ] u 


+ [-cos 0 O + (V afXf0 /g)(a 12 - w 12 )]w + [<V afXf0 /g)(a 13 - w 13 )]q 
+ t U o cos 6 o + W o Sin 9 o + (V a,x,o /g)(a i4 ' W 14 ) ] 0 

+ H^x^ 5 + (V a,x,o/9)(a 15 - w^)]^ + [ (V a> Xf 0 /g) U, 6 - w 16 )]w w 

+ [ (V a,x,o/ g)(b 11 ) ] 6 e + [ (V a f x,o'' g)(b 12 ) l«th 


where a^ ^ to a-j ^ and b-j ^ and 2 represent the constant coefficient terms 
associated with the perturbed-airplane inertial acceleration u in the positive lon- 
gitudinal direction of the body axes, and w^ ^ to 6 represent the constant coef- 
ficient terms associated with the perturbed-wind acceleration u w in the positive 
longitudinal direction of the body axes, derived in appendix B. The term u Q is the 
nominal body longitudinal inertial velocity component, 0 Q is the nominal pitch 
attitude, and u w Q is the nominal wind-velocity component^in the positive longitu- 
dinal direction of the body axes. The terms ^ and V 0 „ _ are the nominal 

a, x, o a, x, o 

longitudinal airspeed and longitudinal airspeed-rate components and g is the accel- 
eration due to gravity. The term u Q is the ^nominal body longitudinal inertial 
acceleration component of the airplane, and u w Q is the nominal wind-shear acceler- 
ation in the positive longitudinal direction of' the body axes and is also derived in 
appendix B. 


The energy-rate signal in equation (C 21 ) was then passed through a second-order 
low-pass filter with appropriate time constants (ref. 11) in order to model the 
energy-rate sensor dynamics. This filter was modeled as follows: 






(C22 ) 


where X^ and are the state and the cutoff frequency associated with the acous- 

tic filter, and X t and a)^. are the state and the cutoff frequency associated with 
the altitude-rate transducer. (From ref. 11, u> f = 1.13 rad/sec and 
= 2.40 rad/sec.) 
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APPENDIX D 


DESIGN MODEL DEFINITION 

This appendix defines the matrix elements of the design model given by equa- 
tions (11) and (12). These equations are repeated here for convenience: 


X 


A 

B 

0 


X 


0 


B w 

• 

u 


0 

0 

0 


o 

+ 

I 

V + 

0 

• 

z 


_ A z 

B z 

0 


z 


0 


0 


1 

** 

I 


1 

o 

o 

\ 



Y 

*0 

= 

0 % 0 


D 

1 


0 0 H z 


Z 


The matrices given in these equations are defined as follows: 


a n 

a 1 2 

a 1 3 

a 1 4 

-a 11 

_a 1 2 

0 

0 

K e ( b i2> 

0 

a 21 

a 22 

a 23 

a 24 

" a 21 

-a 22 

0 

0 

K e (b 22 ) 

0 

a 31 

a 32 

a 33 

a 34 

-a 31 

-a 32 

0 

0 

K e (b 32 ) 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

W 11 

W 1 2 

W 1 3 

W 1 4 

W 1 5 

W 1 6 

0 

0 

0 

0 

W 21 

w 22 

W 23 

W 24 

w 25 

w 26 

0 

0 

0 

0 

k 1 “f 

k 2 a) f 

k 3“f 

k 4“f 

k 5 l0f 

k 6 0)f 

-Wf 

0 

K e o) f ( P 2 ) 

0 

0 

0 

0 

0 

0 

0 

0>t 

-0) t 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

"“e 

0 

h 1 

h 2 

0 

h 4 

0 

0 

0 

0 

0 

0 


(D1 ) 
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a 


z 


21 


0 


H x = 


0 

0 



0 1 0 
0 0 1 
0 0 0 


J 25 


0 



0 0 
0 0 


0 0 0 0 0 
0 0 0 0 0 


0 

0 

0 

1 


h 


1 


0 


h 2 0h 4 0 0 0000 

0 00 0 0 0100 


H 


o ~ 


0 


0 

1 


H 


z 



0 

1 


where 


(D6) 


(D7) 


(D8) 


a z 

Z 21 

- 

^, 0 ^ / V a,o 

(D9) 

a z 

z 22 

(w 0 - 

w w,o^/ V a,o 

(D10) 

3z 25 

* z 21 


(D1 1 ) 

az 26 

“ a ~ 

z 22 


(D1 2) 

V = 

[(u 

- u ) 2 + (w - w ) 2 ] 1/2 

(D13) 

a , o 

L o 

w,o o w,o J 



The K e term is an engine constant that relates engine thrust to throttle position, 
and the value used for the TSRV B-737 was 596 lb/deg. The w 1 j and w 2 j terms in 
equation (D 1 ) are the coefficient terms defined in equations (B22) and (B23), respec- 
tively. The kj and Pj terms in equations ( D 1 ) and (D2) are the coefficient terms 
in equation (Cl 9) of appendix C relating the plant states and controls, respectively, 
to total energy rate. The hj terms in equations (D1 ) and (D6) are the coefficient 
terms for the altitude rate given in equation (Cl 5). The V a Q term is the nominal 
total airspeed. The terms a. - are linearized aerodynamic terms for the airplane. 
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SYMBOLS 


continuous plant-state matrix 
continuous PIF-state matrix 
continuous integrator-state matrix 

2 ^ 2 *] ,A 22 feedforward gain matrices 
continuous plant-control matrix 
continuous PIF-control matrix 
wind-gust disturbance matrix 
PIF wind-gust disturbance matrix 
integrator-state transmission matrix 

coefficient for airspeed and ground-speed feedback gains 

control observation matrix for command model 

error vector used in incremental flight equations 

control matrix for incremental flight equations 

discrete plant-state transition matrix 

discrete PIF-state transition matrix 

discrete closed-loop transition matrix 

discrete command-model state transition matrix 

discrete plant-control matrix 

discrete PIF-control matrix 

discrete command-model control matrix 

discrete PIF disturbance matrix 

acceleration due to gravity, 32.174 ft/sec 

= h Q + h 

PIF-state observation matrix 
command-model state observation matrix 

, plant-state, control-state, and integrator-state observation-matrix 
components comprising the column partitions of the total PIF-state 
observation matrix 


iO> (O' 


column partitions of H corresponding to X, U, and Z states 


VV H Z 

h 

eg 


e / f 


e, t 


ke 


m 


Pe 


I 

i 

J 


K 


K 


e 


k 


altitude of airplane, positive up, ft 
desired-altitude-rate command-model input, ft/sec 
altitude of airplane center of gravity, ft 
specific total energy rate of the airplane, ft/sec 

specific total energy-rate state associated with the pneumatic filter of 
the energy-rate sensor, ft/sec 

specific total energy-rate state associated with the rate transducer of the 
energy-rate sensor, ft/sec 

specific kinetic energy rate, ft/sec 

command-model altitude, positive up, ft 

specific potential energy rate, ft/sec 

identity matrix 

imaginary number (v^T) 

continuous-time cost function 

optimal-feedback gain matrix 

engine constant relating engine thrust to throttle position, 596 lb/deg 
sample integer 


M PIF discrete-cost-function cross-weighting matrix between states and 

controls 

m mass 

N number of samples in discrete performance index 

N PIF output-noise vector 

P A discrete feedforward matrix 

P XZ' P UZ' P ZZ optimal-output feedback solution matrix partitions corresponding to 
X, U, and Z 

Q aircraft total pitch rate, q Q + q, rad/sec (deg/sec in figures) 

continuous PIF-state-weighting matrix 
discrete-cost-function PIF-state-weighting matrix 
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s>> 


*K' Q 0 ,Q Z subsets of Q relating to plant states X, control states U, or 
integrator states Z 

discrete-cost-function PIF control-weighting matrix 
continuous PIF control-weighting matrix 
t time, sec 


U = u Q + u 

U control position-command state vector 

U* ideal control star trajectory 


m 


w 




u 


u w 

u w,E 




Vo 

U W i to 

V 


command-model control vector 


= Vo + "w 

horizontal wind-shear gradient, (ft/sec) /ft 

longitudinal body-axis inertial aircraft-velocity component, positive 
forward, ft/sec 

nominal longitudinal body-axis inertial aircraft-velocity component, 
positive forward, ft/sec 

wind-velocity component in the positive longitudinal direction of the body 
axes, ft/sec 

wind-velocity component in the positive longitudinal direction of the 
Earth-fixed inertial frame, ft/sec 

wind-gust component in the positive longitudinal direction of the body 
axes, ft/sec 

nominal wind-velocity component in the positive longitudinal direction of 
the body axes, ft/sec 

nominal wind-shear acceleration component in the positive longitudinal 
direction of the body axes, ft/sec 

small low-pass cutoff frequency for the horizontal-body-axis wind-velocity 
component, rad/sec 

control rate-command vector 


V a airspeed, ft/sec 

V a,c desired-airspeed command-model input, ft/sec 

v a,m command-model airspeed, ft/sec 

V a,x longitudinal airspeed component in the positive forward direction, ft/sec 
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longitudinal airspeed-rate component in the positive forward direction, 
ft/sec 2 

ground speed, ft/sec 


ground-speed lower limit for the airspeed-ground-speed implementation, 
ft/sec 

measurement-noise covariance matrix 
control-input process-noise covariance matrix 


wind-velocity vector 


= w Q + w 

wind-gust and integrator-noise vector 

total discrete plant-process-noise covariance matrix 

plant-process-noise covariance matrix for random wind-gust disturbances 

= w TT ^ + w TT 
w, o w 

vertical wind-shear gradient, (ft/sec) /ft 

vertical-body-axis inertial aircraft-velocity component, positive down, 
ft/sec 

nominal aircraft inertial velocity component in the positive vertical 
direction of the body axes, ft/sec 

wind-velocity component in the positive vertical direction of the body 
axes, ft/sec 

wind-velocity component in the positive vertical direction of the Earth- 
fixed inertial frame, ft/sec 

wind-gust component in the positive vertical direction of the body axes, 
ft/sec 

nominal wind-gust component in the positive vertical direction of the body 
axes, ft/sec 

nominal wind-acceleration component in the positive vertical direction of 
the body axes, ft/sec 2 

small low-pass cutoff frequency for the vertical-body-axis wind-velocity 
component, rad/sec 

plant-state vector 

PIP-state vector 


ideal state star trajectory 
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X. 

1 


*m 


th 


eg 


*m 


u 


z 

Z E 

ot 

Y 

A 

AT 


6 . . 
ID 


°th 

0 

% 


pneumatic filter state 

initial-condition error covariance matrix 
command-model state vector 
rate transducer state 
throttle-dynamics filter state 

along-track (longitudinal) position of aircraft center of gravity, ft 

initial-condition error vector 

PIF output vector 

command-model output vector 

control-state output vector 

plant-state output vector 

integrator-state output vector 

integrator-state vector 

aircraft altitude, positive down (i.e*, -H) , ft 

aerodynamic angle of attack, deg 

flight-path angle, deg 

perturbation, deviation, or error 

discrete sampling period, sec 

elevator control, deg 

Kronecker delta 

throttle control, deg 

aircraft total pitch attitude, 0 O + 0, rad (deg for plots) 

low-pass-filter cutoff frequency associated with the first-order engine lag 
dynamic model, rad/sec 

low-pass-filter cutoff frequency for the pneumatic filter associated with 
the energy-rate sensor, rad/sec 

low-pass-filter cutoff frequency for the rate transducer associated with 
the energy-rate sensor, rad/sec 
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Subscripts: 


b 

body-axis reference frame 

E 

Earth-fixed reference frame 

k 

present time increment 

k-1 

last time increment 

k+1 

next time increment 

o 

nominal 

T 

total 

Superscripts : 

T 

transpose 

-1 

inverse 

* 

ideal star trajectory 

Special 

notation: 

E { } 

expectation operator 

+DD 

increasing downdraft 

+HW 

increasing headwind 

+TW 

increasing tailwind 

tUD 

increasing updraft 

Abbreviations : 

c.g. 

center of gravity 

DR 

damping ratio 

GM 

gain margin 

MAC 

mean aerodynamic chord 

PM 

phase margin 


A dot over a symbol indicates a first derivative with respect to time, 
dots over a symbol indicate a second derivative with respect to time. 


Double 
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Aircraft parameters Trim values 


Weight, lb 80 000 

c.g., percent MAC . • • . . 25 

Flaps, deg 40 

Glide slope, deg 3 

V_, knots 125 

a 







TABLE II.- CONTINUOUS WEIGHTING AND UNCERTAINTY VALUES 


Variable 

Trim 

States: 


u r ft/sec 

21 1 

w, ft/sec 

5 

q, rad/sec 

0 

0, rad 

-.03 

ft/sec 

0 

w w , ft/sec 

0 

h Sff , ft/sec 

-1 1 

h e t' ft /sec 

-1 1 

x th' deg 

18.5 

h, ft 

1 300 

6 e , deg 

2.6 

6 th' deg 

18.5 

/ h, ft-sec 

0 

/ V , ft 

0 

a 


Controls : 


6 e / deg/sec 

0 

6 th , deg/sec 

0 

Wind: 


U w,g' ft/sec 

0 

w w, g' ft / sec 

0 

Outputs : 


V a , ft/sec 

21 1 

q, rad/sec 

0 

0, rad 

-.03 

h, ft 

1 300 

h, ft/sec 

-1 1 

t , ft/sec 

-11 

<5 e I deg 

2.6 

$th' deg 

18.5 

/ h, ft-sec 

0 

/ V , ft 

0 

a 



Maximum 

perturbation 


Weighting 


0.0004 

.01 

2 500 
10 000 

.0004 
.001 1 
.0044 
.0044 
.04 
.25 
.04 
.04 
.01 
.007 


Uncertainty 


2.5 

.00035 

.004 

3 

1 

1 

.01 

.01 

.01 

.01 


5 
























TABLE III.- FEEDBACK GAINS FOR THE DESIGN WITH ENERGY-RATE FEEDBACK 


Feedback gain for output 


Inputs 


1.73 507 701 3.37 4.39 0.315 -4.87 

-2.00 -9.52 -35.5 -.441 -.397 -.191 .0828 


0.987 0.488 

-1 .52 -.0451 


-0.0093 

-.250 


TABLE IV.- OPEN-LOOP PHUGOID AND SHORT-PERIOD POLES 

Based on SRI International wind-shear-gradient estimations 

for Kennedy airport 


Wind shear 

Phugoid pole 

Short-period pole 

U z = 0.25 (ft/sec)/ft, +HW; 

0.0140 ± i .271 ; 

-0.591 ± il .17; 

W z = 0.47 (ft/sec)/ft, fUD 

DR = 0.0517 

DR = 0.450 

U 2 = 0.39 (ft/sec)/ft, fTW; 

-0.129 ± i .0843; 

-0.690 ± i .909; 

W z = 0.33 (f t/sec)/f t, +DD 

DR = 0.8 38 

DR = 0.604 

U z = 0.39 (ft/sec) /ft, +TW; 

0.104 ± i . 1 08 ; 

-0.620 ± il .06; 

W z = 0.19 (ft/sec) /ft, fUD 

DR = 0.694 

DR = 0.505 
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TABLE VI.- FEEDBACK GAINS FOR THE REDESIGN WITHOUT ENERGY-RATE FEEDBACK 



Input 



Feedback gain for output - 



1.92 510.0 672.0 3.42 

-2.02 -11.6 -19.8 -.45 



/ h 

' V a 

0.478 

-.041 

-0.045 

-.230 



TABLE VII.- ANALYSIS RESULTS OF REDESIGNED CONTROL SYSTEM 
WITHOUT ENERGY-RATE FEEDBACK 


Eigenvalue 

Frequency response for 
Bode plot with - 

Frequency response of 
singular value plot for - 

Closed loop 

DR 

6 e open 

«th open 

Additive 

disturbances 

Multiplicative 

disturbances 


Design case (U_ =0; W„ = 0) 
Z z 


0.617 ± i .531 I I GM, -12.35 dB GM, -15.85 dB 


1 .460 ± i2 .01 0 
-2.400 

1 .860 ± i • 340 
-1 .130 
-.292 

.202 ± i . 278 



PM, 39.50° 


PM, 41 .60 c 


GM / 7 * 28 dB 

GM / 4,45 dB 

( -3.90 dB 

(-9.61 dB 

PM, ±32.97° 

PM, ±39.10° 


Off-design case (U = 0.39 (ft/sec)/ft, +TW; W- = 0.19 (ft/sec)/ft, +UD) 


-0.374 ± i .447 

0.642 

GM / 12.21 dB 

GM, -1 1 .38 dB 

cm/ 5 * 16 dB 

gm/ 3 * 44 dB 

-1 .580 ± i2 .020 

.617 

( -12.39 dB 

PM, 43.30° 

( -3.21 dB 

( -5.77 dB 

-2.430 


PM, 40.20° 


PM, ±25.87° 

PM, ±28.10° 

-2.400 






-1 .130 






-.958 ± i.861 

.744 





-.200 ± i .205 

.698 





-.1 21 






0 






.020 







48 
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Figure 2.- Command -model inputs 






gure 3.- Concluded 




Figure 4.- Concluded 






rformance of desi 




Figure 7.- Performance of energy-rate design with airspeed and ground-speed control 














airspeed and ground-speed control 





1. Report No. 2. Government Accession No. 

NASA TP- 241 2 

3. Recipient's Catalog No. 

4. Title and Subtitle 

TOTAL ENERGY-RATE FEEDBACK FOR AUTOMATIC GLIDE-SLOPE 

5. Report Date 

May 1985 

TRACKING DURING WIND-SHEAR PENETRATION 

6. Performing Organization Code 

505-45-33-04 

7. Author(s) 

Christine M. Belcastro and Aaron J. Ostroff 

8. Performing Organization Report No. 

L- 1 5845 


10. Work Unit No. 

9. Performing Organization Name and Addresa- 


NASA Langley Research Center 
Hampton , VA 23665 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Technical Paper 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


16. Abstract 


Low-altitude wind shear is recognized as an infrequent but significant hazard to all 
aircraft during the take-off and landing phases of flight. A total energy-rate sen- 
sor, which is potentially applicable to this problem, has been developed for measur- 
ing the specific total energy rate of an airplane with respect to the air mass. This 
paper presents control-system designs, both with and without energy-rate feedback, 
for the approach to landing of a transport airplane through a severe-wind-shear and 
gust environment in order to evaluate this application of tne sensor. A system model 
is developed which incorporates wind-shear-dynamics equations with the airplane equa- 
tions of motion, thus allowing the control systems to be analyzed under various wind- 
shear conditions. The control systems are designed using optimal-output feedback and 
are analyzed using frequency-domain control-theory techniques. Control-system per- 
formance is evaluated using a complete nonlinear simulation of the airplane combined 
with a severe-wind-shear and gust data package. This evaluation is concerned with 
control system stability and regulation capability only - operational performance 
aspects are not considered. The analysis and simulation results indicate very 
similar stability and performance characteristics for the two designs. 


17. Key Words (Suggested by Author(s)) 18. Distribution Statement 

Wind shear Unclassified - Unlimited 

Optimal control 
Output feedback 
Energy- rate feedback 

Subject Category 08 

19. Security a assif. (of this report) 20. Security Classif. (of this page) 21. No. of Pages 22. Price 

Unclassified Unclassified 72 A04 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1985 




















